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Abstract 

Quasiparticle approach to dynamics of dark solitons is applied to the case of 
ring solitons. It is shown that the energy conservation law provides the effective 
equations of motion of ring dark solitons for general form of the nonlinear term in 
the generalized nonlinear Schrodinger or Gross-Pitaevskii equation. Analytical 
theory is illustrated by examples of dynamics of ring solitons in light beams 
propagating through a photorefractive medium and in non-uniform condensates 
confined in axially symmetric traps. Analytical results agree very well with the 
results of our numerical simulations. 
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1. Introduction 

As is known, dark solitons are fundamental excitations of nonlinear me- 
dia with modulationally stable background (see review articles [H [2] and refer- 
ences therein). Usually their dynamics is modeled by a generalized nonlinear 
Schrodinger (NLS) or Gross-Pitaevskii (GP) equation with defocusing (repul- 
sive) nonlinearity 

iA + l&f, _ /(|Vf)V = U(r)i>, (1) 

where tf> — ^(r, t) is the field variable which interpretation depends of the phys- 
ical system under consideration. For example, in case of atomic Bose-Einstein 
condensate (BEC) it represents the order parameter ("condensate wave func- 
tion") and then p = \ip\ 2 has a meaning of density of atoms, ip = ^fp exp(i^), 
u = V</> is the condensate's flow velocity, and V(y) is the external potential of 
forces acting on atoms (e.g., the trap potential). In nonlinear optics, the fields 
variable ip has a meaning of the light field strength, |V>| 2 its intensity, t plays the 
role of the coordinate along the light beam, r is the radius-vector in transverse 
directions, and U(r) can be related with non-uniformity of a refractive index. 



Email addresses: kamchOisan.troitsk.ru (A.M. Kamchatnov), 
svyatoslav.korneev@gmail.com (S.V. Korneev) 



Preprint submitted to Elsevier 



October .5, 2010 



In the most common case the function 



/(IVI 2 ) = M 2 (2) 

corresponds to a repulsive interatomic interaction in BEC case or a defocusing 
Kerr nonlinearity in nonlinear optics case. If ip depends on time and one space 
coordinate x only and the medium is uniform (U(r) = 0), then Eq. ([I]) with 
nonlinear term ^ reduces to standard ID NLS equation 

iilH + \^xx ~ = (3) 

with well-known dark soliton solution 

ip = i) s {x - Vt) = | V>o - V 2 tanh ^ Po -V 2 (x~ Vt) +iv}e- ipot , (4) 

where po is the background density and V the soliton's velocity. Naturally, the 
parameters po and V are constant here. 

In many situations we can approximate locally the soliton solution of the 
general equation ([I]) the soliton solution analogous to formula (4) but now with 
slow dependence of po ancl V on space and time coordinates. For example, if a 
dark soliton propagates in BEC confined in a cigar-shaped trap, then its dynam- 
ics can be described in certain approximation by the following generalization of 
Eq. (§: 

# t + \i> xx - M 2 V> = u(x)i>, (5) 

where U(x) = u>qX 2 /2 in case of a harmonic axial potential. As was found in [3], 
such a soliton oscillates in a harmonic trap with the frequency ui = oj /\/2, and 
it was shown in [3] that this result can be obtained most simply if one considers 
a dark soliton as a quasiparticle with the energy 

£{ V) = l(p -V^ (6) 

which expression follows easily from the Hamiltonian form of Eq. |3]). Then 
conservation of energy 

e(V) = const (7) 
combined with the Thomas-Fermi (TF) distribution of the background density 

Po(x) = jji — U(x), p = const, (8) 

yields at once the equation for V = dX/dt, where X(t) is a coordinate of 
"position" of the soliton: 



dJC- 2 
~dt 



+ U{x) = const. (9) 



This approach was developed in more detail in Ref. [5] where, in particular, its 
generality was illustrated by application to the case of nonlinearity f(p) = p 2 . 
General formulae for the case of arbitrary function f(p) (provided that the 
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condition of modulation stability of the background is satisfied) were derived 
later in Ref. [§]. It is clear that this approach is correct as long as width ~ 
1 / \/po — V 2 01 the soliton is much less than a characteristic length at which the 
background parameters change considerably so that a "collective" coordinate 
X(t) is a well-defined variable. 

This approach admits also the form of Hamilton equations which describe 
dynamics of the soliton quasiparticle. Indeed, as was shown in [7] (see also [5]), 
the momentum of the soliton Q is equal to 



p = -2V^ Po - V 2 + 2p arctan(v/po - V 2 /V) (10) 
and the equation 

dp dt K ' 

holds which evidently is one of the Hamilton equations for the quasiparticle with 
the Hamiltonian (JsJ) . In case of the dark soliton moving through non- uniform 
BEC confined in a trap created by the potential U(x) the density po is a function 
of x related with U(x) by Eq. Hence the Hamiltonian becomes a function 
of the position X of the soliton, e(X, V) = e(po(X), V), where po(X) is related 
with the potential U(X) by the equation ^ with x replaced by the coordinate 



X of the soliton's position. It is clear that the equation (11) does not change 
provided the derivative in the left-hand side is taken for fixed X. Then the 
second Hamilton equation 

dp _ fde{X,V) 



can be readily reduced to 



it = -\r^l (12) 



dV__dU(X) 

1 dt - ox [U) 

that is to the Newton equation with the effective mass of the dark soliton quasi- 
particle equal to the doubled atomic mass as it was first discovered in [3] (see 
also 01 El). 

In this Letter we show that the approach of Refs [H El E] can be easily 
generalized to multi-dimensional situations of ring or spherical solitons. For 
example, if a ring dark soliton evolves in a nonlinear medium, then the speed of 
its radius R(t), V = dR/dt, depends not only on time t, but also on R. If the 
width of the soliton is much less than its radius, then the conservation law of 
the total energy in the first approximation reads 

2irRe(p (R), dR/dt) = const. (14) 

Thus, we arrive at a differential equation for R(t) which can be easily solved 
in various situations. We shall consider here by this method several typical 
examples. Naturally, in the simplest case of Kerr nonlinearity and uniform 
background our approach reproduces well-known results of Ref. [8j. In actual 
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optics experiments (see, e.g., [9]) on ring solitons, the photorefractive media with 
f(p) = p/(l + 7p) were often used, and in the next section we shall consider 
evolution of ring dark solitons in such a medium. Evolution of ring dark solitons 
in BEC confined in axially symmetric traps was considered in Refs. [TOl E] in 
various approximations. We show in section 3 that our approach gives a very 
good approximation (different from ones discussed earlier) which agrees quite 
well with direct numerical solution of the GP equation up to the moment of 
decay of ring solitons to vortices due to their snake instability. 



2. Dynamics of dark ring solitons in a uniform background 

First, we shall consider the case of uniform background which density (in- 
tensity) is constant, po — const. Although new results will be obtained for 
a non-Kcrr nonlincarity function /(/?), we shall begin with the Kerr-like case 
f(p) = p which illustrates very clearly the method. 

2.1. Kerr-like nonlinearity 

In case of Kerr-like nonlinearity (|2|), i.e. f(p) = p, the energy of a dark 
soliton per unit of its length is given by Eq. Hence the total energy of a 
ring dark soliton with radius R is equal to 2irRs and its conservation law yields 
the equation 

R(p -R 2 ^ 2 ^R (p -V 2 ^ 2 , (15) 

where R = dR/dt and Rq and Vq are the initial radius and velocity of the dark 
soliton, correspondingly. This equation was derived first in a different form in 
[3] with the use of perturbation theory. Its differentiation with respect to time 
t gives the Newton-like equation 

R=^(p -R 2 ) = R °\~ V °\ (16) 
3R H ' 3i? 5 / 3 v ' 

which has been recently obtained in [12 as a particular case of equations of 



contour dynamics of dark solitons. Variable R can be excluded from ( 16 1 with 

(17) 



the help of Eq. ( 15 ) so that one obtains 

dV (po - V 2 f 



dt 3i?o(po-V 2 )3/ 2 ' 

and elementary integration of this equation yields implicit dependence of V on 
time 

V(3p - 2V 2 ) = V (3p - 2V 2 ) - t/R 

(p -y2 )3 /2 (p - 1/2)3/2 ■ I > 

This equation can be solved with respect to V 2 in a closed form (see [T^]) and 
this analytical solution agrees very well with the results of direct numerical 



solution of the GP equation. Thus, interpretation of Eq. ( 15 ) as the energy 



conservation law permits us to reproduce very simply the known results. 
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2.2. General form of nonlinearity 

As was shown in [5] , in general case the energy of a dark soliton per its unit 
length can be expressed as 

irow^ (19) 

where 



fPO 

Qo(p)=8p [f(po)-f(p')W, (20) 
J P 

Q{p) = Qv{p)-W 2 { Po - P )\ (21) 

and the minimal value p m in the distribution of the density p across the dark 
soliton dip is related with its velocity V by the equation 

V 2 = Q v(p m ) (22) 
4(po~P m ) 2 

Hence, e can be considered as a function of po and V 2 and, as above, conservation 
of the total energy of the ring soliton reads 

Re(p ,R 2 ) = R e( Po ,V 2 )- (23) 

This is a differential equation for the function R — R(t) with given initial 
conditions R = Rq, V = Vq at t = 0, and its solution describes the evolution 
of the ring dark soliton. To obtain more definite results, one should specify the 
nonlinearity function f(p). 

2.3. Evolution of ring dark solitons in media with photorefractive nonlinearity 

In nonlinear optics experiments [9] the nonlinear effects can be modeled by 
the function 

f(p) = — (24) 
1 + 1P 

which takes into account saturation of the nonlinear response, i.e. f{p) —5-1/7 = 
const as p — > 00. In this case 



Qo(p) = ^(lni^° + ^^). 
Y \ 1 + 7/3 l+7/Oo 



Substitution of this expression into Eqs. (19l and (22| yields the function e(po, V ) 



and the resulting differential equation ( 23 ) can be solved numerically. Example 



of such a solution is shown in Fig. 1 for the choice of parameters: po = 1, 7 = 1, 
i?o = 30, Vq = 0.2. The solid line represents the results of our approximate 
theory and they are compared with the exact numerical solution of the general- 
ized NLS equation |l]) with / given by Eq. (24 1 and U = 0; values of the dark 



soliton radius for different moments of time are shown by crosses. As we see, 
agreement is quite good. 
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Figure 1: Dependence of the radius of the ring dark soliton on time for the case of photorefrac- 
tive nonlinearity: solid line corresponds to analytical theory and crosses to numerical solution 
of the generalized NLS equation. 



3. Evolution of ring dark solitons in a non-uniform background 

As a typical example, we shall consider the BEC confined in an axially 
symmetric trap. As above, we assume that cylindrical symmetry persists during 
evolution so that dynamics of a ring dark soliton reduces to varying of its radius 
R with time t. However, now the background density po is not uniform and 
depends on the radial coordinate r. In the Thomas-Fermi (TF) approximation 
and for quadratic nonlinearity f(p) — p this dependence is given by 

Po (r) =n- U(r) (26) 

and in the case of a harmonic trap U(r) = WqT 2 /2 this dependence can be 
written as 

p (r) = ^ 2 (l 2 -r% (27) 

where I is the radius of the condensate's distribution related with the chemical 
potential p and the number of particles N (per unit axial length) by the formulae 

P=\^l\ N = ^u$l i . (28) 

Therefore the energy conservation law now reads 

Re( Po (R),R 2 ) = R e(p Q (R ),V 2 )- (29) 
We shall study it in some detail for a particular case of a harmonic trap when p$ 



is given by Eq. (27) and of two-particle repulsive interatomic interaction when 



6 



f(p) = p and e(V) is given by Eq. (J6J). Then Eq. (29) takes the form 



R 



3/2 



i?o(po-^ 2 ) 3/2 , (30) 



where po = uj 2 (1 2 — Rq)/2 is the density at the initial radius R$ of the soliton. 
This equation can be rewritten as 

dR\ 2 _ 1 2 2 2 R 2 Q /3 (p -V 2 ) 

~^ o{l ~ R) W* ' (31) 

Its differentiation with respect to time t gives the Newton-like equation 
d 2 R ^ Rl / \p -Vj) 

W~~Y R + 3R^ ■ (32) 

The first term in the right-hand side corresponds to oscillations of the soliton 
radius with the frequency i.e. it describes the effects of non-uniformity 

of the condensate discovered in [3]. The second term in the right-hand side 



of (32 1 describes the effects of the curvature of the soliton and coincides with 



Eq. (16). In our adiabatic approximation both effects are weak and therefore 



they are added in the combined equation (32 1. Integration of this equation 
reduces to a quadrature 

"°'= ± r,. f ..... °= ^''t; v °\ m 



V2 J Ro /l Vl-W 2 - aW - 2 /3 ' wg/8/3 

which actually solves the problem. In Fig. 2 we have compared this analytical 
solution with the exact numerical solution of Eq. |l]) for the choice of the pa- 
rameters N = 3486.5, fi = 1, I = 50, u = 0.03, R = 10.0, V = 0.229. At 
the moment of time t c = 79.3 the ring soliton reaches its turning point with 
the maximal value of the radius R ma x = 37.25; after this moment its radius 
decreases and, as it was already indicated in Ref. [TU], the snake instability 
starts to develop, so that the dark soliton evolves into necklace of vortices. This 
process is illustrated in Fig. 3. In Fig. 2 the radius of the necklace is shown by 
circles and formal analytical solution corresponds here to a dashed line. As a 
measure of deviation of the structure from the ring soliton we have introduced 
the quantity 



5 P =^P 2 -{-p)\ (34) 

where the averages p and p 2 are taken along the circle of the minimal values of 
the density in the structure. Its dependence on time is shown in Fig. 4. It can 
be fitted by the exponential growth of bp: 

<5pw 8.198 • 10" 9 exp(0.123t) (35) 

The coefficient can be represented as 8.198 • 10~ 9 = 0.0065 • exp(-0.123 • 110) 
which corresponds to the numerically introduced disturbance about ~ 0.01 
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Figure 2: Dependence of the radius of the ring dark soliton on time in case of its propagation 
along TF distribution of density BEC confined in a harmonic trap: solid line corresponds to 
analytical theory and crosses to numerical solution of the GP equation. 



which starts to grow up after the moment of time about t s» 110, that is af- 
ter the turning point what can be explained as a result of very slow change 
of parameters of the soliton in vicinity of the turning point so that the most 
unstable mode has enough time for its evolution without strong competition 
with other modes. The growth rate constant T — 0.123 can be compared with 
that calculated according to the theory [13] of snake instability of plane dark 
solitons for given depth and background density. If we take these parameters 
for the moment t = 110, when the soliton radius equals to R = 30.9, then we 
get the maximal growth rate T = 0.09 for the unstable mode with the wave 
number k = 0.39 in qualitative agreement with the numerically found value. It 
is worth noting that we have here the number of wavelengths along the ring of 
the soliton equal to Rk = 12 in agreement with the number vortices formed at 
the nonlinear stage of development of the instability. Thus, the snake instabil- 
ity theory agrees at least qualitatively with the process observed in numerical 
simulations. 

4. Conclusion 

The idea that bright solitons behave as "particles" goes back to the begin- 
ning of modern nonlinear physics. However, it was not so obvious that it could 
be applied to dark solitons since their behavior depends also on the background 
evolution. Konotop and Pitaevskii noticed 0] that if background is a stationary 
solution of the GP equation in the TF approximation, then dynamics of dark 
solitons is determined by conservation of their energy. We have shown in this 
Letter that this idea can be extended to dynamics of ring solitons. This simple 
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Figure 3: Transformation of a ring dark soliton into necklace of vortices: density plots at the 
moments of time (a) f = 120 with a slightly disturbed dark soliton and (b) t = 130 with well 
developed necklace of vortices. 



approach can be easily generalized to arbitrary form of the nonlinear repulsive 
self-action what is important for applications to nonlinear optics. Examples dis- 
cussed in this Letter show that this approach gives accurate enough description 
of motion of ring dark solitons up to the moment of their decay to vortices due 
to snake instability. 
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